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Low-temperature thermodynamics of the classical frustrated ferromagnetic spin chain near the 
ferromagnet-helimagnet transition point is studied by means of mapping to the continuum limit. The 
calculation of the partition function and spin correlation function is reduced to quantum problem 
of a particle in potential well. It is shown that exactly at the transition point the correlation length 

f^ , behaves as T~^" and the magnetic susceptibility diverges as T"*'-^ in the low-temperature limit. 

^Nj ' Corresponding numerical factors for the correlation length and the susceptibility is calculated. It 

is shown that the low-temperature susceptibility in the helical phase near the transition point has 
a maximum at some temperature. Such behavior as well as the location and the magnitude of the 
maximum as a function of deviation from the transition point are in agreement with that observed 

r~^ ' in several materials described by the quantum s = 1/2 version of this model. 

(N' 
o ; i. introduction 

^ . 

Lately, there has been considerable interest in low-dimensional spin models that exhibit frustration [l|. A very 
_j interesting class of such systems with unique physical properties is chain compounds consisting of edge-sharing CuO^ 
.^^ units [2tI|- The frustration in these compounds arises from the competition of the ferromagnetic (F) interaction 
. ' J I of nearest neighbor (NN) spins and the antiferromagnetic (AF) next-nearest-neighbor (NNN) interaction J2. An 
"trt ! appropriate model describing the magnetic properties of such copper oxides is so called F-AF spin chain model, the 
rj Hamiltonian of which has a form 

'^ . H = Ji 2_^ Sn • S„+i + J2 2_^ S„ • S„+2 (1) 

O , where Ji < and J2 > 0. 

, ^, This model is characterized by the frustration parameter a = J2 / 1 Ji I . The ground state properties of the quantum 

s = 1/2 F-AF chain have been intensively studied last years |8l-[l6|. It is known that the ground state of model ([T]) 

y—^ ■ is ferromagnetic for a < 1/4. At a = 1/4 the quantum phase transition to the incommensurate singlet phase with 

^ ' helical spin correlations takes place. Remarkably, this transition point does not depend on a spin value, including the 

'O . classical limit s -> 00. 

However, the influence of the frustration on low-temperature thermodynamics is less studied, especially in the 

vicinity of the ferromagnet-helimagnet transition point. It is of a particular importance to study this problem, 

because edge-sharing cuprates with a ~ 1/4 (for example, Li2CuZrOi, Rb2Cu2Mo-iO 12) are of special interest \u\ . 

t^ ' Unfortunately, at present the low-temperature thermodynamics of quantum s = 1/2 model ([1]) at a ^ can be 

studied only either by using of numerical calculations of finite chains or by approximate methods. On the other hand, 

the classical version of model ([Ij can be studied by analytical methods giving exact results at T — > 0. Of course, 

• • , the question arises about the relation of these results (in particular, for the susceptibility) to those of the quantum 

. !^ ' niodel. It is known [18l420l | that quantum and classical ferromagnetic chains [a — 0) have universal low-temperature 

K> , behavior. As was noted in Ref.[l9| the physical reason of this universality is the consequence of the fact that the 

5^ ' correlation length at T ^^ is larger than de Broglie wavelength of the spin waves. This property is inherent in the 

Cu ^ frustrated ferromagnet too. Though such universality for the frustrated ferromagnetic chains is not strictly checked 

at present, one can expect that the universality holds on for the F-AF chain as well. Therefore, the study of classical 

model ll]) can be useful for the understanding of the low-temperature properties of the quantum F-AF chains. 

At zero temperature classical model ([T|) has long range-order (LRO) for all values of a: the ferromagnetic LRO 
at a < 1/4 and the helical one at a > 1/4. At finite temperature the LRO is destroyed by thermal fluctuations 
and thermodynamic quantities have singular behavior at T — >■ 0. In particular, the zero-field magnetic susceptibility 
X diverges. For the ID Heisenberg ferromagnet (HF) x — 2|Ji|/3T^ ^21,]. At < a < 1/4 the susceptibihty is 
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X = 2(1 - 4a) |Ji| /ST^. This behavior of x is similar to that for the quantum s = 1/2 F-AF model ^^. The value 
XT^ vanishes at the transition point. As it was noted in Ref.[22| this fact indicates the change in the critical exponent. 

In this paper we focus on the low-temperature behavior of the classical F-AF chain near the ferromagnet-helimagnet 
transition point. At first we consider the case a = 1/4, i.e. the F-AF model exactly at the transition point. This 
problem is interesting on its own account, because the spectrum of low-energy excitations is proportional to fc^ rather 
than fc^ as for the HF model. It means that the critical exponents characterizing the low-temperature behavior of 
thermodynamic quantities at a = 1/4 can be different from those for the HF chain. Besides, the method developed 
for the study of the transition point can be generalized to investigate the vicinity of the transition point. 

At the ferromagnet-helimagnet transition point a = 1/4 it is convenient to rewrite Hamiltonian ([T]) in the form 
[231: 

i^=^^(S„+i-2S„ + S„_i)2 (2) 

In Eq.([2]) we put | Ji| = 1 and omit unessential constant. 

In the classical approximation the spin operators S„ are replaced by the classical vectors Sn of the unit length. 
In what follows we use the continuum approach replacing Sn by the classical vector field s{x) with slowly varying 
orientations, so that 

Sn+l — '^Sn + Sn-1 ~ TT^ (3) 

where the lattice constant is chosen as unit length. 

In the low-temperature limit the thermal fluctuations are weak, so that neighbor spins are directed almost parallel 
and continuum approach ^ is justified. Using the continuum approximation, Hamiltonian ^ goes over into the 
energy functional of the vector field s[x): 

1 r . fd^Y 



This energy functional is a starting point of the investigations of model ([T]) at a = 1/4. The partition function is a 
functional integral over all configurations of the vector field on a ring of length L 



Z^jDsi.)e.^^--J^ dx^-J| (5) 

It is useful to scale the spatial variable as 

C = 2T^^^x (6) 

Then, the partition function takes the dimensionless form 

Z./z,.(Oexp{-/^«(0)] (7) 

where the rescaled system length is A = 2T^^^L. The partition function ([7]) and the correlation function (s(Z) ■ s(0)) 
are the objects of the current study. 

The paper is organized as follows. In Sec. II we consider the planar version of spin model Q at a = 1/4. For 
this more simple model we demonstrate the technique of the calculation of the correlation function. We show that 
the thermodynamics of this classical one-dimensional model reduces to the zero-dimensional quantum problem of a 
particle in a potential well. In Sec. HI the classical continuum F-AF model at the transition point is studied. In 
this case the partition function describes a quantum particle in an axially symmetrical potential well. We obtain the 
exact expressions for the susceptibility and the structure factor. In Sec. IV the behavior of the uniform susceptibility 
in the helical phase at a ^ 1/4 is studied and compared with the experimental data for the edge-shared compounds 
and with the results for the quantum s — 1/2 model. The conclusions are summarized in Sec.V. and the Appendix 
contains some technical aspects of the calculation of the correlation function. 



II. PLANAR SPIN CASE AT a = 1/4 
A. Partition function 

We begin our investigation of the thermodynamics in the transition point with a more simple planar spin version 
of model ([5]), when all spin vectors lie in one plane and have only two components: 

s(e) = (sin0(O,cos0(O) (8) 

Such order of study is methodically justified, because the technique of the correlation function calculation is similar 
for both planar and original three-component spin models, but it is easier to demonstrate on the simple planar model. 
In terms of 6'(^) the Hamiltonian transforms to 

JA' f<PeV fdo^^ 

and the partition function becomes 



dej \de) ' Kd^j ^^^ 



Z^ I D9{C) exp J - / dC{ 



(10) 



where the prime denotes the space derivatives d/d^. 

In general, when one deals with the field theory containing the second or higher order derivatives one has to 
follow the Ostrogradski prescription [24]. However, as will be demonstrated below, in our case we can avoid such 
complications and calculate the partition function and correlation functions in a more simple way. 

Since the Hamiltonian contains only derivatives of the field 0{£^), the partition function can be rewritten in terms 
of a new field 



Z = J Dqi^e^pl-J rfe 



{<l" + l')} (11) 

J [^ Jo 

where 



,(0 ^ « (12) 

d^ 

To calculate the partition function we utilize well-known equivalence of the n-dimensional statistical field theory 
with the (n — l)-dimensional quantum field theory. It is obvious in advance that partition function dill) describes a 
quantum particle in a potential well U{q) = q^ at 'temperature' 1/A. However, we will follow all intermediate steps, 
because we will need them in the subsequent calculations of the correlation function. 

The transition amplitude (or propagator) of a particle located initially at g(0) = g^, and finally at q{t) = qf takes 
the form of a path integral 

{qf\e-''" \q,) ^ p Dq{t)exp[i 1^ dtL{q,q)\ (13) 

Then, imposing the periodic boundary conditions qf — qi — q and integrating over g, we obtain the partition 
function in a form 

itH 



Z^ j dq{q\e-^'"\q) (14) 

In our case we replace ^ by an imaginary time ^ ^ it and partition function (jlip takes the form of a path integral 
of a quantum particle in a potential well: 



Z= f Dq{t)expli f dtLo{q,q)\ 



(15) 
J \^ Jo J 

where the Lagrangian is 

Lo^q^- q" (16) 
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FIG. 1: Wave function tpo{q) for the fowest level £o = 0.4208 of planar spin model. 

The monicntuni p is p ^ 2q and the Hamiltonian is 



Ho = ^P^ + q^ 
The corresponding Schrodinger equation describes a quantum anharmonic oscillator: 

-4d? + '^^ = '^ 

The spectrum of equation (|18p is calculated numerically: 

Ea = 0.4208; 1.508; 2. 96... 

and the wave function ipo{q) for the lowest level Eq = 0.4208 is shown in Fig.l. 
Now the exponent of Hq can be represented as follows: 



,-AH( 



"-Eiv^")' 



-ASa 



(V'a 



(17) 



(18) 



(19) 



(20) 



and the partition function becomes 



Zcx / dq{q\e-^"\q 



)-E' 



-A£a 



(21) 



As expected, we obtain the partition function of a quantum anharmonic oscillator at 'temperature' 1/A. In the 
thermodynamic limit A = {2T^/^L) — > cx) only the lowest eigenvalue Eq = 0.4208 gives contribution to the partition 
function, 



Z ^e 



-Aeo 



(22) 



B. Correlation function 



The scalar product of vector fields located on distance I can be written as 



?(0 • s(0) = cos [61(0 - 6'(0)] = 



')'{x)dx 



= » 



exp ( i / q{£,)d£, 



(23) 



where ^ = 2T'^/H. 

Then, the correlation function can be represented as a ratio of two functional integrals: 



(s{l)-m) = ^nZc] (24) 

where denominator Z is already calculated (see Eq. (f2T|) ) and 

Zc= f Dqex^pl- f di (g'2 + q^) + z f qd^ I (25) 

The latter path integral is clearly divided on two parts [0, fi] and [fi, A] and Zc can be represented as 

Zc= dqodqfj,Zi{qo,q^)Z2iqi,,qo) (26) 

where the propagators Zi and Z2 are 

Zi((7o,<Zm) = P ^'i''^^\ j '^ dtLi{q,q)\ (27) 

Z2{q^,qx) ^ f^Dqexpi-f dtLo{q,q)] (28) 

and periodic boundary condition qx = qo is applied. 
The Lagrangian Lq is given by Ea. ([T5)) and 

Li=f-q'^ + iq (29) 

The propagator Z2 is calculated straightforward using Eqs. ([T^ and (PH)) : 

Z2 = ^e-(^-A')^» (gj^^) (^^|go) (30) 

But the propagator Zi requires special treatment, because Li and the corresponding quantum Hamiltonian Hi are 
non-Hermitian: 

^.^-ii^ + .^-.. (31) 

Non-Hermitian operator Hi can be represented as 

7?! = ^ 77q |Uq) {Va\ (32) 

a 

where |uq,) and \va) are eigenfunctions of direct and conjugate eigenvalue equations: 

Hi \Ua) = r/Q |Ua) 

H\ K) = < |z;„) (33) 

The normalization conditions for \ua) and \va) are 

(t;„|w^) = {Ua\vp) = baS (34) 

Equations (I33p for Hamiltonian (I3ip transform to each other by complex conjugation operation. This implies that 
eigenfunctions in Ea.p3p satisfy the relation \va) = |m*). Thus, we need to solve only one of the differential equations 
([55)1 . The numerical calculations show that all eigenvalues of Eg. (155)1 are real and positive. A few lowest eigenvalues 
are presented in Eq.([35l). 

7;„ = 0.6472; 1.517; 2.99... (35) 
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FIG. 2: Real and imaginary parts of the wave function uo{q) for the lowest level rjo = 0.6472 of planar spin model. 



Real and imaginary parts of uo{q) for the lowest level rjo — 0.6472 are shown in Fig. 2. 
Now, the propagator Zi can be expressed through the solutions of Eg. (1551) . 



Zi = (go I e-^"^ \q^) = ^e-^"" {qo\u^) {v^ 



?M/ 



(36) 



where we used the identity 



-f^H 



^=Ek)' 



-fJ-Vc 



(37) 



Then, substituting Egs. ipn]) and (pS)) into Eg. ip^ and integrating over qo,qf_i we obtain 



^c = E' 



-A'('?o-e/3)„-Ae/3 



(V'/3|Ma> 



«„8 



The main contribution in the thermodynamic limit A ^ oo is given by the lowest value eg and Zc reduces to 



^c^e-^^«Ee""^''°"'°^('/'o|wa)' 



(38) 



(39) 



Substituting /i = 2T^'^l into Eq. (p9|) . we find from Eq. (f24| the correlation function 



(40) 



The correlation length is governed by the lowest eigenvalue rjQ and equals to 

1 



In 



2ri/3(,yo - eo) 



= 2.2T"^/^ 



(41) 



So, the low-temperature behavior of the correlation length is different from the HF model, where Ic ^ T ^. 
Now, the structure factor can be also calculated 



S{k) = 2K / dle'^^ {s{l) ■ s(0)) = 2^ (V'olua)^ 
Jo 



2rV3(^^_eo) 
4r2/3(r;„ - eo)2 + fc2 



(42) 



In the low-temperature limit, the expansion of the structure factor for any k ^ T^/^ has the form 

^(^) = -^- E(^" - ^O) (^Ol"")' - T^ E (^0\Uaf {Va ~ S^f + . . . (43) 

a. a 

The first term in Eq. (|43l) is zero, because 

E(^" - So) (^o|ua)' - (V^ol ^1 - -ffo l^o) = -« (V'ol g l^o) = (44) 

and IfAo) is even function of q. 

Therefore, the structure factor is given by the second term in Eq.(j43|). which can be calculated exactly 

E(^" - so)^ (^ol^a)' = (V'ol (^1 - %f + \ [[^0, i^i] , i^i] IV^o) ^-\ (45) 

Therefore, the low-temperature asymptotic of the structure factor is 



4T 
5(fc) = ^ (46) 

Hence, the susceptibility x(fc) foi' ^ 3> T^''^ remains finite in the low-temperature limit: 

X(*) ^ fi 4 (47) 

The fact that x(^) ~ ^~^ (instead of fc~^ for HF chain) is a consequence of the fact that the excitation spectrum 
becomes ~ fc** at the transition point. 

For fc = the structure factor (|^ diverges at T ^ as 

5(0)-r-i/3yMffi)! (48) 

The sum in Ea. (|48|) is calculated numerically and gives the factor w 5.36. Therefore, the magnetic susceptibility 
behaves as T"**/^ 

X(0)^aa = g| (49) 

In conclusion of this section we emphasize that the exact calculation of the correlation function for the planar spin 
model demonstrates that the critical exponents at the transition point of the F-AF chain can differ from that for the 
HF chain. 

III. CLASSICAL SPIN MODEL AT THE TRANSITION POINT 
A. Partition function 

The calculation of the correlation functions for the classical three-component spin model ([5]) is to a large extent 
similar to the planar spin case and to avoid duplications we will often refer to the previous section. So, in this section 
we consider the continuous model described by energy functional ^ where three-component vector field s(^) has the 
constraint s^(C) = 1- 

Since in the partition function the integration occurs over all possible spin configurations, we are free to choose 
any local coordinate system. It is convenient to choose it so that the Z axis at the point ^ is directed along the spin 
vector s(f ), so that the spin vector s(f ) — (0,0,1). 

Let us introduce a new vector field 

ds 
q{0^ -J7^{qx,qy,qz) (50) 



The constraint ^{^) — 1 converts to the relations for q{^): 



4 - < (51) 



where the prime denotes the space derivatives d/d^. 
Then, the Hamilton function in Eq.([7]) transforms to 



0y = (|)'='z;^+'?;^ + ('^^+9^)^ (52) 

Here we see that the constraint ^ — 1 effectively eliminates the Z component of q from the Hamilton function. 
Therefore, henceforth we deal with the Qx and qy components of the vector field qonly, and we denote a two-component 
vector field by q(f) = {qx,qy)- 

The partition function in terms of q(^) takes the form: 

Z = Joqexpl-J^ de(q" + q')| (53) 

Similar to the planar spin case, we treat the partition function as path integral (jlSp for the quantum mechanics of 
a single particle with the Hamiltonian 

^o = -iA + q4 (54) 

where A = 9^ + 9Ms two-dimensional Laplace operator. Hamiltonian ((54|) commutes with the z-componcnt of the 
angular momentum l^ and eigenstates tp{q) of the corresponding Schrodinger equation 

Hoijj = Sip (55) 

are divided to subspaces of azimuthal quantum numbers Iz — 0, ±1, ±2 . . .. 

Thus, the wave function ipi^ (q) describes a particle with the azimuthal quantum number Iz in 2D axially symmetrical 
potential well U{q) = q'^. Numerical solution of Eg. ([55]) gives the lowest levels for Iz = 0, ±1, ±2 

e^{l^ = 0) = 0.9305; 3. 78; 7.44... 
ea{lz = ±1) = 2. 14; 5.48; 9.44... 
ea{lz = ±2) = 3. 54; 7.27; 11. 5... (56) 

The wave function V'o('?) for the lowest eigenvalue eo = 0.9305 is shown in Fig. 3. 

Now, by the analogy with the planar spin case Eas. ([20l) - (|21l) we obtain the partition function as a sum of exponents 
over all quantum numbers a and Iz'- 



J^e-'^'-^^ (57) 



In the thermodynamic limit A ^- cxa only the lowest level Sq = 0.9305 survives and the partition function is 

Z = e"^^" (58) 

B. Correlation function 

Since we work with the local coordinate system directed so that the spin vector is directed along the Z axis, in 
order to find the scalar product of spin vectors s{l) ■ s(0) we need to express the vector s(0) in the local coordinate 
system located at the point x = I. This can be represented as a chain of successive rotations describing the trajectory 
s{x) [see the Appendix]: 

sx=i{0) ^explt f Qdx] Sx=QiO) (59) 
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FIG. 3: Wave functions tp{q) (solid line) for the lowest level £o = 0.9305 and U3{q) (dashed line) for the level 770 = 1.4113. 



where fl — [a ■ w(a;)) can be expressed through the vector q{x) as 

„ {(f-[qxcf]) Qxq'y-qyq'^ 

il = = (7yq.j. - ffxqy + 0-2 



cp 



ll+ly 



(60) 



In this equation we used relations (j5ip for the Z component of q and (f . 

Then, taking into account that the spin vectors are directed along the Z axis of local coordinate system Sx=o{0) 
Sx=i{l) — (0, 0, 1), the scalar product of spin vectors s{l) ■ s(0) becomes: 



s{l) • s(0) = ( 1 ) cxp I i /" r2(q, q')dx j I 



(61) 



Similar to the planar case, the correlation function can be written as a ratio of two path integrals ([24)) . where the 
denominator Z is given by Eq. ()58p and the numerator represents the following path integral: 



JDqcxpl-J d^ (q'2 + q") + t J"^ r!(q, q')d^ I 



(62) 



Then, repeating the arguments presented in Eas. (j25l) - pOI) . we arrive at the problem of the calculation of the 
propagator: 



9m 



i^qexp t / dtLiiq, q) \ = {qo\ e"^^^ \q^) 

qo i Jo ) 

where the Lagrangian and the corresponding quantum Hamiltonian are 



(63) 



ii = q - q - iayqx + icJxqy -cTz , , , 

Ix + % 

1 A i ^z ^ '^'^z^z 

Hi = --A + q H —2 \- layqx - la^qy 

Substituting aa from Eqs. (|A.5p . the Hamiltonian takes the matrix form 



Hi 



^0 + -^ 2^ilz Qx 



2q 



yHz Ho + j^ qy 



-qx 



4q2 
Qy Ho 



(64) 
(65) 

(66) 
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where Hq is defined by Ea. ((54)) . 

Operator Hi is non-Hermitian and the exponent of Hi can be represented as 



e-pffi 



^|w„)e-''''°(z;„| (67) 



where three-component eigenvectors u{qx,qy) = (""1,^2, "3) and v{qx,qy) = (wi,W2,W3) satisfy the corresponding 
eigenvalue equations 

Hiu — r/u 

hIv = ri*v (68) 



The normalization conditions are 



(Valup) = {Ua\v/3) = (5a,/3 (69) 



Making the same procedure for non-Hermitian operators as for the planar spin case we obtain the correlation 
function in a form 

(s(0 • m) = 3ff E (V'ol^a^a) (^a.alV'o) e-^^"""^") (70) 

a 

Only the eigenfunctions Jms^q) and {v^^a) are present in the above equation, because according to Eg. dFTI) we need 
only the element (3, 3) of the resultant matrix. 

Since the wave function ipQ has zero angular momentum, then only the sector Iz — of Eqs. (l68p gives the contribution 
to the correlation function. In this sector the wave functions depend only on q — |q|, the Hamiltonian is simplified so 
that we have to solve a pair (instead of three) differential equations for u^^q) and (j){q) — (qxUi + qyU2)/q'- 



l(P<j) 1 d<t) 1 
4 dq^ Aq dq Aq 

1 d'^u^ 1 du3 4 



7 3::j-z;:7r + z;2'^ + '?'^ + 9"3 - r,<j, 



4 dq^ Aq dq 



q U3-q(l) = 77U3 (71) 



The conjugate eigenvalue problem for v in Eg. (1551) transforms to exactly the same differential equations ([7T|) for 
the functions v^{q) and x*{<l) = ^(fe^i +qyV2)/q- Therefore, the function V3{q) is found from the solution of Eq. ([7T|) 
by the relation v^lq) = u^{q) and the the normalization conditions (I69|) transform to 

(w3,al"3,/3) - (CI-^/s) = Sa,f3 (72) 

One can see that equations ((7T|) describe a two-level system in an axially symmetric potential well U{q) — q'^, where 
two levels with angular momenta Iz = and Iz = 1 are coupled by non-Hermitian transition operator. The spectrum 
of this system of equations turns out to be real and positive as for the planar case and a few lowest levels are: 

7/„ = 1.4113; 1.83; 3.98... (73) 

The reality of the spectrum t^q, has an important consequence: the correlation function (|70l) decays on large distances 
without oscillations. 

The wave function u-i{q) for the lowest level 770 = 1.4113 is shown in Fig. 3. As follows from Fig. 3 the behavior of 
the functions ipo{q) and 1*3(9) are similar. Therefore, the main contribution to the correlation function, the structure 
factor and the susceptibility is given by the level 770- 

The correlation function and the structure factor are given by equations (j40p . (1421) with the substitution u^^a for 
Ua and eigenvalues presented in Egs. dSB]) . ([75]). Therefore, the correlation length defined by the lowest eigenvalue r]o 
behaves similar to the planar spin case r-^ T"^/'^, but the numerical factor is different: 

Ic = 1.04T^^/^ (74) 

The low-temperature asymptotic of static structure factor S{k) and the susceptibility x(fc) for k ^ T^/^ is calculated 
in a similar way as for the planar spin case (see Eg. (l43p ) resulting in 

S{k) = f 

X(fc) = 3^ (75) 



11 



For k = the structure factor S{0) is defined by Ea. p8|) with the corresponding eigenfunctions and eigenvalues. 
The numerical calculation of the sum gives 



^(0)=T-V3^(V^lK^^|^ (76) 



Thus, we have arrived at the final result for the magnetic susceptibility of the classical spin model at the transition 
point: 

,., S{0) 1.07 
We see that the planar spin model gives correct critical exponent for the magnetic susceptibility. 

IV. LOW-TEMPERATURE SUSCEPTIBILITY IN THE HELICAL PHASE 

In the preceding sections we have considered the low-temperature thermodynamics at the transition point a — 1/4. 
Likewise it is possible to extend the developed method for the study of the vicinity of the transition point. In this 
case an additional term 2{a — l/4)s'^ appears in energy functional ^, which after rescaling of spatial variable ([5]) 
forms the scaling variable 

Especially interesting is to study the vicinity of the transition point when values {a — 1/4) and T are small but 
the parameter 7 is fixed. All the steps in derivation of the expressions for the correlation function are exactly the 
same as was done for the transition point and only the form of the potential well in the corresponding differential 
equations ((7T|) is modified: U{q) = q^ ~ A'yq^. Numerical solution of the corresponding differential equations allows 
to find the correlation function, the static structure factor S{k) and the susceptibility as a function of 7. As follows 
from Eqs. (|76|) and (|77|) the susceptibility can be rewritten as 

/(7)^^!Z(2L (79) 

^ T4/3 (a -1/4)2 V»; 

Thus, the normalized susceptibility X — x(Q! ~ 1/4)^ is a universal function of 7. 

On the ferromagnetic side of the transition point (a < 1/4) the susceptibility diverges at T — ^ 0, but the exponent 
changes from 4/3 to 2, so that the susceptibility becomes x ^ (1/4 — a)/T'^. 

The behavior of the susceptibility in the helical phase is more interesting. The dependence of the uniform suscep- 
tibility X on the normalized temperature x = 7"^/^ _ T/{a — 1/4)^^'^ in the helical phase is shown in Fig. 4. The 
characteristic features of this dependence are the maximum oi x ^^ ^ — ^m and the finite value of x at T — >■ 0. The 
latter fact is a classical effect and can be destroyed by quantum fluctuations. In the quantum s = 1/2 F-AF model 
the ground state is believed to be gapped (though the gap can be extremely small [Ef) and so x — ^ at T -> 0. In 
real materials interchain interactions cause the 3D spiral LRO and the susceptibility can remain finite at T = 0. 

The obtained dependence x(^ is in a qualitative agreement with those observed in the edge-shared compounds 
with a close to 1/4 {Li2Cu02 [3, Rb2Cu2Mo30i2 Q, Li2CuZr04 [l3)- As follows from Eq.^ the location of 
maximum of xiT) is at Tm ^ {a— l/4)'^/2 ^nd Xm '^ (a — 1/4)"^^ i.e. with the increase of a the maximum of x shifts 
to higher temperatures and the magnitude of the maximum Xm decreases. These dependencies of Tm and Xm on the 
frustration parameter a are also in accord with the experimental observations [17] . 

The dependence xi^) agrees also with the numerical results for the uniform susceptibility obtained by TMRG and 
exact diagonalization methods for the quantum F-AF chain with s = 1/2 [l^jllll- The dependence Tm{a) is similar 
to that obtained by the TMRG calculations in Ref.[l4|. Thus, the classical model catches the physics of quantum 
spin systems in the helical phase and, therefore, the developed method for the classical spin model represents a useful 
tool for the investigation of the low-temperature behavior of the quantum systems. 

V. SUMMARY 

We have obtained the exact results for the low-temperature thermodynamics of the classical F-AF model at the 
frustration parameter a = 1/4, where the ground state phase transition from the ferromagnetic to the helical phase 
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FIG. 4: Normalized susceptibility x(a ^ 1/4)^ as a function of normalized temperature T/(q — 1/4) 
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occurs. The main result relates to the behavior of the zero-field susceptibility x and the correlation length Ic- It is 
shown that the critical exponents of x and Ic are changed from 2 to 4/3 and from 1 to 1/3 correspondingly, when 
a — >■ 1/4 from the ferromagnetic side. We note that the critical exponents are the same both for the classical spin 
model and for the planar model. In the present paper we have considered a continuum version of the model ^. 
However, the low-temperature asymptotes of x and Ic of the continuum and lattice models coincide. In fact, lattice 
model ([U can be studied on a base of the transfer matrix method adapted in Ref.[25] to include the NNN interaction. 
We have shown [26] that the corresponding transfer-integral equations are reduced at T — !■ to the differential 
equations (1551) . (j7ip with the same eigenfunctions and eigenvalues, so that the exact low-temperature asymptotic of 
the susceptibility coincides with Eg. ((77)) . Besides, the structure factor iS'(fc) for the lattice model is given by Eq. p2)) 
under substitution of 2(1 — cos/c) for k^ in the dominator of Eq. (|42l) . At the same time, calculations in a frame of 
path integral method are essentially simpler and clear in comparison with those in the transfer-integral approach. 

It is interesting to compare the exact expression for the susceptibility at the transition point with the results found 
by approximate methods. One of this methods is the modified spin-wave theory (MSWT) proposed by Takahashi [27| 
to extend the spin- wave theory to the low-dimensional spin systems without LRO. Another approximate approach is 
the expansion of the thermodynamics of the n- vector classical model in powers in 1/n [2^ [2^ (usually, up to the first 
order). Remarkably, both methods give the true critical exponent 4/3 for the susceptibility x = cT~^/^. However, the 
numerical coefficient c differs from the exact one. For example, the MSWT result is c = 1.19. The 1/n expansion for 
the classical spin model ^ gives c = 0.560 in the zeroth order and c — 0.897 in the first order in 1/n. The comparison 
of these values of c with the exact one shows that these approximate methods give a satisfactory agreement (within 
10 — 20%) with the exact coefficient. It is worth to note that the MSWT gives the exact low-temperature asymptotic 
of X = 2/3T^ for the classical ferromagnetic chain {a — 0) [21]. Moreover, the MSWT gives the exact result for x at 
T — > for the quantum ferromagnetic chain with s — 1/2 as well. As was noted in Introduction, the quantum and 
the classical ferromagnetic chains have universal low-temperature properties. The low-temperature susceptibility of 
the ferromagnetic chain is described by the scaling function universal for any value of spin s. For the F-AF model at 
a = 1/4 there is no rigorous proof of such universality, though MSWT confirms this hypothesis. If this universality 
holds in the case a = 1/4 then we expect that the quantum F-AF chain has the same critical exponent of x as in the 
classical model. 

We have also considered the low-temperature thermodynamics in the vicinity of the transition point. In this case 
the properties of the system are governed by the scaling parameter ^ — {a — l/4)/r^/'^. On the ferromagnetic side 
of the transition point {a < 1/4) the susceptibility transforms smoothly to x ^ (1/4 — ct)/T'^ at 7 ^ — cxd, which 
describes the change in the exponent at T ^' 0. The susceptibility in the helical phase for a fixed value of a ^ 1/4 
behaves as x ^ (« — l/4)~^ at 7 ^- cx), which means that the uniform susceptibility remains finite at T — > 0. Besides, 
we found that the susceptibility in the helical phase has a maximum at some temperature T^ ~ (a — l/4)~^'^. The 
presence of maximum of the dependence x(r) as well as the location and the magnitude of this maximum as a function 
of the deviation from the transition point (a — 1/4) are in agreement with that observed in several materials described 
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by the quantum s = 1/2 version of this model and with the numerical results for the s = 1/2 model. 
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Appendix 



Let us consider a particle moving along a trajectory f{t). At any moment the particle motion can be represented 
as an instant rotation around some instantaneous axis of rotation with definite angular velocity. The radius of the 
instant rotation p is expressed through the centripetal part fc of the acceleration vector r: 



r 



(A.l) 



The angular velocity is 



UJ = — = 
P 



[f X r\\ 
r 



(A.2) 



where we took into account that (t^, ■ f) ~ 0. 

Since the rotation takes place in the (r , r) plane, the instantaneous axis of rotation is directed along [r x r\ . Therefore, 
if we define the 'local' coordinate system associated with the moving particle so that the instant coordinate axes are 
directed along the vectors f, fc and [r x f], then the instant change in the local coordinate system is expressed by the 
rotation matrix 



R = 



UJ 



exp [i [a ■ uj) 



where 
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(A.3) 
(A.4) 

(A.5) 



This implies that if a given fixed vector s has components sto in the local coordinates corresponding to time to, 
then the components of this vector in the local coordinates corresponding to time ii can be represented as a chain of 
successive rotations: 



S(i = exp 



((? • uj{t)) dt 



to 



Stn 



(A.6) 
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